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ABSTRACT: The objective of this study was to estimate genetic parameters that are to be used for across-herd genetic evaluations of 
seed stock pigs at GGP level. Performance data with pedigree information collected from swine breeder farms in Korea were provided 
by Korea Animal Improvement Association (AlAK). Performance data were composed of final body weights at test days and ultrasound 
measures of back fat thickness (BF), rib eye area (EMA) and retail cut percentage (RCP). Breeds of swine tested were Landrace, 
Yorkshire and Duroc. Days to 90 kg body weight (DAYS90) were estimated with linear function of age and ADG calculated from body 
weights at test days. Ultrasound measures were taken with A-mode ultrasound scanners by trained technicians. Number of performance 
records after censoring outliers and keeping records pigs only bom from year 2000 were of 78,068 Duroc pigs, 101,821 Landrace pigs 
and 281,421 Yorkshire pigs. Models included contemporary groups defined by the same herd and the same seasons of births of the same 
year, which was regarded as fixed along with the effect of sex for all traits and body weight at test day as a hnear covariate for 
ultrasound measures. REML estimation was processed with REMLF90 program. HeritabUity estimates were 0.40, 0.32, 0.21 0.39 for 
DAYS90, ADG, BF, EMA, RCP, respectively for Duroc population. Respective heritability estimates for Landrace population were 
0.43, 0.41, 0.22, and 0.43 and for Yorkshire population were 0.36, 0.38, 0.22, and 0.42. Genetic correlation coefficients of DAYS90 with 
BF, EMA, or RCP were estimated to be 0.00 to 0.09, -0.15 to -0.25, 0.22 to 0.28, respectively for three breeds populations. Genetic 
correlation coefficients estimated between BF and EMA was -0.33 to -0.39. Genetic correlation coefficient estimated between BF and 
RCP was high and negative (-0.78 to -0.85) but the environmental correlation coefficients between these two traits was medium and 
negative (near -0.35), which describes a highly correlated genetic response to selection on one or the other of these traits. Genetic Trends 
of all three breeds tend to be towards bigger EMA or greater RCP and shorter DAYS90 especially from generations bom after year 2000. 
(Key Words: Swine, Duroc, Landrace, Yorkshire, Heritability, Growth, Ultrasound) 



INTRODUCTION 

The most popular pig breeds in Korea are Landrace and 
Yorkshire for maternal lines and Duroc for paternal line. 
There are some, but not in large scale as Duroc, minor 
paternal line breeds used in small scale, which are Berkshire 
or Hampshire. Recent industrial needs for genetic 
evaluation of pigs for selection are increasing with global 
trends for seed stock pig patient. In Korean system of pig 
breeding, grand-grandparent (GGP) stockers maintain 
purebred animals. Grandparent stockers make cross 
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between Landrace and Yorkshire lines (PS) for distribution 
to commercial farmers. Conamercial feedlot operators 
usually adopt paternal lines as terminal sires to make three- 
way crossbreds. 

The objective of this study was to estimate genetic and 
environmental variances of days to 90 kg body weight and 
ultrasound meat traits and covariances between these in 
purebred pig populations. 

MATERIAL AND METHODS 

Performance testing 

Pigs weaned at around three to four weeks of ages, 
when they weigh around 35 kg, were grouped together by 
sexes into performance testing until they reach 70 to 110 kg 
body weight. At least four pigs per litter were selected for 
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testing. At the time of testing, body weight was measiu'ed 
and ultrasound scanning on the back was made with A- 
mode (amplitude mode) scanners. Back fat thicknesses on 
the shoulder (on the 4th thoracic vertebrae), mid-back (on 
the last thoracic vertebrae) and loin (on the last lumbar 
vertebrae) were averaged. Days to 90 kg body weight 
(DAYS90) was calculated according to Korean Swine 
Performance Recording Standards (KSPRS), which 
assumes body weight at birth as 1 kg. 

Days to 90 kg = ((90 kg-body weight at the test) 
x(age at the test in days-38)) 
/(body weight at the test) 

Eye muscle area (EMA) was estimated with A-mode 
ultrasound measures scanned at 5 cm ventral to the dorsal 
point of the last thoracic vertebrae. Average of three backfat 
thickness measurement was adjusted for body weight by the 
formula suggested by KSPRS. 

BF = scanned average backfat+((90 kg-body weight at 
the test)xscanned average backfat)/(body weight at the test- 
11.34) 

EMA from ultrasound scanner were used directly for 
analyses without body weight adjustments. Retail cut 
percentage (RCP) was estimated also from A-mode scanner 
by the function embedded and used for analyses without 
body weight adjustment as well. 

Models 

Models for parameter estimation were separate for 
growth performance (DAYS90) and ultrasound scanned 
measurements. 

(DAYS90) = nH-CG-Hsex-Hu-ne 

(BF, EMA, RCP) = ^-HCG-Hsex-HPi(finwt)-Hu-He 

Here, DAYS90: the days to 90 kg body weight 
calculated by the formula 

ADG: average daily gain (g) from birth to test day 

BF: average of back fat thicknesses (mm) measured at 



three points 

EMA: eye muscle area (cm^, estimated from A-mode 
ultrasound scanner) 

RCP: retail cut percentage (%, also estimated from A- 
mode ultrasoimd scanner) 

CG: contemporary group effect (year and season of 
birth-farm) 

sex: sex categories (1 = gilts, 2 = boars) 

finwt: final weight (body weight at the end day of 
performance testing, kg) 

And |Li, Pi, u, e are population mean, linear regression 
coefficient, breeding value of the animal (BV), temporary 
environmental effect, respectively. 

EM REML estimates of the genetic and environmental 
variances and covariances between traits were estimated by 
remlf90 program (Mizstal, 2002). For random effects (u and 
e), all the genetic or environmental covariances between 
traits were assumed not zero. But the covariances between 
genetic and environmental effects were assumed to be all 
zeros. 

Standard error estimates of the heritabiUty estimates 
were approximated by the scalar modification of the method 
suggested by Fischer et al. (2004). The derivation of the 
variance of heritability estimate used in our study is as 
follows. 

Var(h^) = Var(sVs^) 

^(sVs'p)'-(Var(s'G)/sVVar(s^)/s%) 
-2Cov(s'g, s'p)/(s'g • s'p)) 

RESULTS AND DISCUSSION 

Table 1 show the data structure for analyses by breeds. 
There were tendency to have greater number of boars for 
testing in Duroc population and greater number of gilts in 
Landrace or Yorkshire populations. This might be the 
farmers' preference to keep animals of paternal or maternal 
use after completion of performance testing. 

Estimates of fixed effects are summarized in Table 2 by 
breed groups. Estimates for sex effects were calculated 
from averaging out the other fixed effects (effects of 



Table 1. Pcdi 


gree and perl'oniiaiice data slruclure by breeds 












Total No. 

records 


No. inbred 

animals 


Average 






No. levels 




Breed 


inbreeding 
coefl'icienl ('/i) 


Animal 




Sex 


Contemporary group 


Duroc 


78,068 


45,290 


4.0588 


86,794 


Boar 
Gilt 


47,230 
30,637 


2,104 


Landrace 


101,821 


62,363 


4.7432 


112,742 


Boar 
Gilt 


24,970 
76,753 


2,316 


Yorkshire 


281,412 


169,634 


3.4281 


302,314 


Boar 
Gilt 


42,332 
238,539 


2,985 
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Table 2. Solutions and least squares 


means for the fixed effects from least squares models 






Breed 


No. records 


Classes 




Estimates for fixed effects 






DAYS90 


BF 


EMA 


RCP 


Landrace 


76,753 


Gilt 


149.24 


13.03 


47.80 


56.28 




24,970 


Boar 


143.18 


11.92 


47.25 


56.23 






y * 




-0.0031 


0.0052 


0.0008 


Yorkshire 


238,539 


Gilt 


146.68 


12.98 


46.88 


55.73 




42,332 


Boar 


139.99 


11.49 


46.11 


55.90 










-0.0045 


0.0114 


-0.0017 


Duroc 


30,637 


Gilt 


147.45 


12.98 


46.30 


54.13 




47,230 


Boar 


143.03 


11.65 


45.67 


55.00 










-0.0042 


-0.0028 


0.0027 



* Partial regression coefficient of trait values on body weight on test day. 



contemporary groups) and applying average end weight 
multiplied by corresponding regression coefficients. Days 
required for gilts to reach 90 kg body weight (DAYS90) 
were larger than for boars by about 6 d in Landrace or in 
Yorkshire breeds and by about 4 d in Duroc breed. Gilts 
tended to have thicker backfat (BF) and larger eye muscle 
area (EMA) but smaller retail cut percentage (RCP) 
estimated by ultrasound scanner on the final day of the test. 
Partial regression coefficients of BF, EMA, and RCP on 
body weight at the test were all close to zero but were all 
significantly different from zero as they were estimated in 
preliminary analyses with least squares generalized linear 
models (data not shown). Estany et al. (2002) observed that, 
in the lines of Landrace pigs, a line selected for bigger litter 
size had thicker subcutaneous backfat than a random mating 
control line. However, no such typical trend was found in 
our study populations by the differences among breed 
groups. Even if direct comparisons between breed groups 
are not possible because of no connectedness between them, 
Duroc pigs had almost similar but a little longer DAYS90 
and a little thicker BF than Yorkshire pigs. Landrace pigs 
tended to have longer DAYS90 and thicker BF than 
Yorkshire or Duroc pigs, which might somehow support the 
evidence of maternal breed selection for bigger litter size as 
was found out of the experiments of Estany et al. (2002). 
However, Petry et al. (2004) could not find any strong 
evidence of correlated responses in growth performance or 
in carcass characteristics as responses to selection for Utter 
size. Therefore, paternal or maternal breed differences 
weekly observed in our study populations cannot be 
understood as differences in the directions of selection by 
breeds. 

Variance components 

Genetic parameter estimates in Table 3 show that days 
to 90 body weight (DAYS90) or pork yield predictor traits, 
BF, EMA and RCP estimated from ultrasound scanners are 
all moderately to highly heritable similar to the values in 



the Uterature (Chen et al., 2002; Noguera et al., 2002; Kim 
et al., 2004; Arrango et al., 2005). Hicks et al. (1999), on 
the other hand, reported much higher heritabihty estimates 
(0.43 for ADG and 0.69 for BF) in Japanese Large White 
population. Suzuki et al. (2006) reported higher heritability 
estimates in Duroc pigs; 0.51 for ADG, 0.48 for EMA and 
0.73 for BF when they measured pigs at around 105 kg 
body weight and scanned with B-mode. Heritabihty 
estimates of Li and Kennedy (1994)'s in Canadian 
populations were lower for days to 100kg body weight 
(0.26 to 0.32) and higher for back fat thickness (0.50 to 
0.55) than the results from our studies. The heritability 
estimates in this study were similar among three breed 
groups. 

BF's in all three breeds were lowly but positively 



Table 3. Genetic variances (diagonals), co variances (above 
diagonals), heritability estimates (h^±SE) and genetic correlation 
coefficients (below diagonals) 



Vg/Tg 


DAYS90 


BF 


EMA 


RCP 


Duroc 










DAYS90 


38.15 


0.5394 


-1.929 


1.993 


BF 


0.089 


0.9632 


-0.7844 


-1.147 


EMA 


-0.149 


-0.382 


4.384 


1.349 


RCP 


0.216 


-0.783 


0.431 


2.23 


h^ 


0.4+0.0000 


0.32+0.0004 


0.21±0.0000 


0.39+0.0002 


Landrace 










DAYS90 


43.54 


0.2873 


-3.712 


2.667 


BF 


0.033 


1.695 


-0.9891 


-1.845 


EMA 


-0.248 


-0.334 


5.168 


1.367 


RCP 


0.239 


-0.837 


0.355 


2.866 


h^ 


0.43+0.0004 


0.41+0.0003 


0.22±0.0002 


0.43+0.0004 


Yorkshire 










DAYS90 


32.31 


2.59E-02 


-2.18 


2.75 


BF 


0.004 


1.603 


-0.9975 


-1.882 


EMA 


-0.17 


-0.349 


5.104 


1.574 


RCP 


0.277 


-0.851 


0.399 


3.052 


h^ 


0.36+0.0000 


0.38+0.0000 


0.22±0.0000 


0.42+0.0000 
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correlated genetically or environmentally with DAYS90. 
However, in the study results of Kim et al. (2004), who 
estimated parameters in a Korean company, a subpopulation 
of our study populations showed negative genetic and 
phenotypic correlation between these two traits in all three 
breeds (Duroc, Landrace and Large White) as similar to the 
estimates by Li and Kennedy (1994) in Canadian 
populations. In the report of Chen et al. (2002), the genetic 
and phenotypic correlations between days to 113.5 kg and 
backfat thickness were estimated all to be negative but close 
to zero with STAGE test data for all breeds. Genetic 
correlation between BF and EMA were all moderate and 
negative, which was in agreement with the estimates by 
Chen et al. (2002). However, the genetic correlation 
between EMA and DAYS90 were all high and negative, 
which again controvert to the estimates by Chen et al. 
(2002) whose estimates were low but positive. This 
contradiction in correlation structure between DAYS90 and 
ultrasound measures of BF or EMA would be the difference 
in linear adjustment of ultrasound, which, in our case, no 
prior adjustment was apphed instead of using body weight 
as a covariate in the estimation models. RCP estimated from 
ultrasound was positively correlated genetically with 
DAYS90 or with EMA and negatively correlated with BF in 
all breeds. 

Environmental correlation coefficients shown in Table 4 
were similar to genetic correlation coefficients (Table 3) in 
direction of association between traits. But the 
environmental correlation coefficients between BF and 
EMA (-0.016 to -0.115) or between BF and RCP (-0.277 to 
-0.363) were estimated to be much smaller in magnitude 
than genetic correlation coefficients between those (-0.334 
to -0.382 between BF and EMA, and -0.783 to -0.851 
between BF and RCP). This is supportive of the idea that 

Table 4. Environmental variances (diagonals), covariances (above 
diagonals), and environmental correlation coefficients (below 
diagonals) 



Ve/Ie 


DAYS90 


BF 


EMA 


RCP 


Duroc 










DAYS90 


58.28 


1.275 


-6.121 


3.817 


BF 


0.116 


2.075 


-0.3591 


-0.745 


EMA 


-0.195 


-0.061 


16.88 


1.728 


RCP 


0.268 


-0.277 


0.225 


3.482 


Landrace 










DAYS90 


56.56 


1.091 


-6.579 


3.858 


BF 


0.092 


2.46 


-0.6255 


-1.068 


EMA 


-0.203 


-0.093 


18.52 


1.789 


RCP 


0.265 


-0.352 


0.215 


3.75 


Yorkshire 










DAYS90 


56.26 


0.5626 


-6.774 


3.782 


BF 


0.046 


2.665 


-0.7925 


-1.213 


EMA 


-0.214 


-0.115 


17.89 


2.184 


RCP 


0.246 


-0.363 


0.252 


4.187 



RCP, that is predicted from ultrasound scanned values of 
backfat thickness and eye muscle are, is mostly determined 
by the genetic potential of BF. And EMA or RCP are more 
associated with genotypes of BF than with feeding 
environments. Genotypes of EMA or environmental effects 
on EMA are positively affecting RCP but as much as those 
ofBF 

Genetic trends 

Genetic trends described as average breeding values 
over the years of births in three major breeds studied are 
illustrated in Figure 1 along with average inbreeding 
coefficients (aveF). Average inbreeding coefficient tended 
to rise rapidly from around 1998. But it became rather 
stable from 2000 to 2011 with a big drop in 2004 especially 
in Duroc herds. This might indicate a possible explanation 
that most Duroc seed stock farmers had replaced boars with 
imported ones once for every five years by the year 2004. 
Among the traits evaluated, it seems that selection of boars 
and gilts of all three breeds had been made towards larger 
EMA. Genetic potential for DAYS90 tended to increase by 
around year of birth 2000, which is an unfavorable direction 
in terms of daily body weight gain or feed efficiency. But 
then after, DAYS90 tend to decrease sharply. Genetic 
potential of Duroc herds for BF tended to decrease and that 
for RCP tended to increase especially the generations of 
pigs born from 2000. However, in Landrace or Yorkshire 
herds, genetic potentials for BF tend not to change, but 
those of RCP tend to increase steadily from generations 
bom from early 1990's. These genetic changes seem to 
follow the similar pattern to the trend in average inbreeding 
coefficients. 

In our study, BF was modeled with linear covariate 
function of final body weight along with the other 
ultrasound scaimed measure traits (EMA and RCP) while 
BF was pre-adjusted for body weight already. But still 
significant amount of variation due to regression term of 
body weight was left. Partial regression of BF on body 
weight was negative for all breed groups. And the genetic 
and environmental correlation coefficients between 
DAYS90 and BF were all positive, different from literature 
estimates. This might indicate that pre-adjustment of 
ultrasound backfat measure for body weight at the time of 
measure was not enough to adjust for differential growth of 
pigs around the test end day. Negative partial regression 
coefficients of BF on body weight mean that pigs weighing 
heavier than average tend to have thinner backfat. This is 
also implied by the positive genetic correlation between 
DAYS90 and BF shown from our study even if the 
correlation coefficients were small near zero. A possible 
explanation to the genetic correlation structure between 
traits is that pigs by the end day of performance test 
weighing around 90 kg were immature to show fat growth. 
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Figure 1. Genetic trends of days to 90 kg and ultrasound traits(scales on left vertical axis) with average inbreeding coefficients (aveF, 
scales on right vertical axis) by the year of births(BF = average EBV of backfat thickness; DAYS90 = average EBV of days to 90 kg 
body weight; EMA = average EBV of eye muscle area; RCP = average EBV of retail cut percentage). 



However, significance of partial regression of BF on body 
weight implies that there are still considerable individual 
variations with regards to differential body growth. 
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